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ABSTRACT 

We present the Galactic Arecibo L-Band Feed Array Hi (GALFA-Hl) survey, and its first full data 
release (DRl). GALFA-Hl is a high resolution (~4'), large area (13000 deg^), high spectral resolution 
(0.18 km s-^), wide band (-700 < ^lsr < +700 km s"^) survey of the Galactic interstellar medium 
in the 21-cm line hyperfine transition of neutral hydrogen conducted at Arecibo Observatory. Typical 
noise levels are 80 mK RMS in an integrated 1 km s~^ channel. GALFA-Hl is a dramatic step forward 
in high-resolution, large-area Galactic Hi surveys, and we compare GALFA-Hl to past, present, and 
future Galactic Hi surveys. We describe in detail new techniques we have developed to reduce these 
data in the presence of fixed pattern noise, gain variation, and inconsistent beam shapes, and we 
show how we have largely mitigated these effects. We present our first full data release, covering 7520 
square degrees of sky and representing 3046 hours of integration time, and discuss the details of these 
data. 

Subject headings: ISM:Hl; atlases; surveys; radio emission lines 



1. INTRODUCTION 

In 2004 the Arecibo L-Band Feed Array (ALFA) 
was installed on the Arecibo 305-meter radio telescope. 
While Arecibo's status as the largest single-dish antenna 
in the world had always made it a premier instrument 
in cm-wave radio astronomy, its ability to map large ar- 
eas of sky had been hindered by the fact that it could 
only observe one position at a time. One can consider 
the original L-band wide receiver, which covers a similar 
frequency range to ALFA, as a 'single-pixel camera' with 
a single ~ 3.5' beam. ALFA, by comparison, monitors 
7 positions at once, thus increasing the mapping speed 
of Arecibo by a factor of 7. This tremendous boost in 
efficiency led many different science teams to consider 
for the first time the possibility of mapping huge swaths 
of sky with Arecibo to create a new generation of high- 
fidelity, high-resolution sky atlases. These groups were 
divided by topic: pulsar (P-ALFA; Cordes 2008), extra- 
galactic (E-ALFA ; Giovanelli et al. 2005, Henning et al. 
;2010[|Auld et al.|i2006, Freudlin^ et al 2008) , extrater- 
restrial intelligence ( [Korpela et al. ' 200l)| , and Galactic 
science (GALEA). The GALl^A consortium itself consists 
of three independent groups, or sub-consortia, divided by 
the radiative processes they measure: radio recombina- 
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■ et al.|20'Q5|, co ntinuum 
& 'i'aylor||2009|) and the 



tion lines (GALFA-RRL ; [Terzian 
radiation (GALFACTS; jGuram fc Taylor||^ 
21 cm line transition of H i (GALI^A-Hl; Goldsmith 
[GALFA-Consortium||2QQ3l this work). 

Successful large Hi surveys of the Galaxy have been 
undertaken in the past, perhaps be ginning in the mod- 



ern er a wit h the pioneering work o f Weaver fc Williams 
(|1973|) and|Heiles fc Habing| ( |197 4[). More recently there 



surveys. 



The 



has been an explosion of Galactic Hi 
Leide n- Argentine-Bonn (LAB) Survey ([Kalberla et al. 



2005) is a full sky, stray-radiation corrected map of the 



Galaxy in Hi at relatively low angular resolutio n (36^ 



The International Galactic Pla ne Survey (IGPS; Taylor 
"et al.^'2003[ |Stil et al.|[2QQ6l and [McClure - Griffiths et al. 
2005) is a compilation of relatively high resolution (^Tj 
interferometric observations of t he Galac tic plane. The 
Galactic Ah Sky Survey ( GAS S ; [M^ChlFe- Griffiths etll] 
2009), is a half-sky intermediate resolution (15') survey 
in the southern hemisphere. GALFA-Hl combines the 
large-area approach of GASS and LAB with resolution 
only somewhat coarser than the IGPS. These qualities, 
along with high spectral resolution and superb fidelity, 
allows GALFA-Hl to contribute to a whole new range of 
science in and around our Galaxy. A number of other sur- 
veys planned or in progress are also tackling the imaging 
of large areas of the Galactic sky in Hi. We have tabu- 
lated the main parameters of these completed, ongoing, 
and planned surveys in Table [l] for ease of comparison. 

The G ALFA-Hl survey was designed to study a broad 
range of topics under the umbrella of the neutral hy- 
drogen. Galactic interstellar medium (ISM). To date, 
many different scientific endeavors have been undertaken 
with GALFA-Hl. We have studied gas in the Galactic 
halo, w ith projects based on observations of high- velocity 
clouds ([Peek et al.||2QQ7l iStanimirovic et aL||2008[ [Peek 



et al. 2009, and Hsu et al. "2010|) as well "as constramts 



on t he scale of gas outflows from stars in globular clus- 



ters 



van Loon et al. 2009). We have also investigated 



2 



gas much clos er to home, with studi es of turbulence in 
the local ISM ( Chepurnov et al.|2QlQ ) a nd studies of cold 
neutral medium clouds close to the sun (Peek et al. 2010). 
Compact clouds, small Hi clouds distinct trom dit tuse 
background emission, are a focus of investigation (Be- 



gum et al. 2010 



Saul et al. 2010). GALFA-Hi has also 



been leveraged to study the complexities of th e Galac 



tic plane in Hi ( iKrco et al.| |2007| |Gibson||2010[ )^_as_w^ 
as the details oi galaxies in the Local Group (Putman 



et al.||2009| |Grcevich et al.||2008| . Excluding the Galac 
tic plane work, all ot these projects are based on data 
that is presented to the public in this work. 

In this paper we present a description of the GALFA- 
Hl survey in progress and the first data release (DRl). 
DRl covers 7520 square degrees of sky, utilizing 3046 
hours of Arecibo data taken with ALFA. A large fraction 
of this data was taken commensally with other projects, 
allowing for the efficient use of telescope time. The sur- 
vey covers the velocity range —700 < '^lsr < +700 km 
s~^ at 0.184 km s~^ spectral resolution. Due to the dif- 
fering integration times over the survey footprint (see 
r igure 1 , RMS noise in a 1 km s ^ channel ranges from 
140 mK to 60 mK, with a median value of 80 mK. The 
final effective beamsize is 3.9 arcminutes in RA and 4.1 
arcminutes in Declination, slightly larger than the raw 
ALFA beamsize. Figure [2] shows an area covering roughly 
50% of the DRl data release, and demonstrates the in- 
credible fidelity and spatial dynamic range of the sur- 
vey. The reduced data are available to the public at 
https://purcell.ssl.berkeley.edu/, and the raw data are 
also available upon request. 
This paper begins by describing the overall signal chain 
and the observing techniques we employed in gather- 
ing the data (J3|. We also describe the way the data were 
processed (§4F7focusing on some of the new techniques 
we developea to process this unique data set, along with 
an example of the effect of the methods we employ in 
a sample region (Q. We discuss the data presented in 
DRl (^, and outline our plans for the future of the 
ongoingGALFA-Hl survey (Q. 

2. THE SYSTEM 

2.1. The AO 305 m 

The GALFA-Hl survey is conducted on the Arecibo 
Observatory 305-meter telescope, south of Arecibo, 
Puerto Rico, the central telescope of the National As- 
tronomy and Ionosphere Center (NAIC). The primary 
reflector can be described as the 305 m diameter cap of 
a 265 m radius sphere. Suspended above the primary 
reflector is a triangular platform that carries a circular 
track upon which the azimuth arm rotates. The azimuth 
arm carries a geodetic dome (also called the Gregorian 
dome), which houses secondary and tertiary optics to 
correct for spherical aberration. The entire system can 
track objects from Declination -1° 20' to 38° 02' with 
~ accuracy. ALFA is located on a rotating turret in- 
side the dome at the focal plane of the telescope. The 
details of the 305 m are discussed in the Arecibo Techni- 
cal and Operations Memos Series (ATOMS) and on the 
NAIC website, http://www.naic.edu. 

2.2. ALFA 



ALFA is a 7-beam feed array designed to allow Arecibo 
to perform surveys at a 7-fold increased rate. The re- 
ceiver frontend was built under contract to the Australia 
Telescope National Facility of the Australian Common- 
wealth Scientific and Industrial Research Organization, 
based upon the successes of the P arkes multi-beam re- 



The beams are ar- 
I a beam in the mid- 



ceiver ( |Staveley-Smith et al. 1996 

ranged on the sky in a hexagon wit 

die (see Figure [3|). Each beam is slightly elliptical along 
the zenith-angle direction, approximately 3.3' x 3.8' full- 
width half-maximum (FWHM) on the sky, with a gain of 
~11 Jy/K for beam and ~8.5 Jy/K for beams 1-6. Each 
beam is divided into independent dual-linear polariza- 
tions, for a total of 14 independent signals. ALFA can be 
tuned from 1225 to 1525 MHz, which includes the hyper- 
fine transition of Hi at 1420.406 MHz. For more design 
and implementation details consult the ALFA technical 
memo series, http://www.naic.edu/alfa/ineinos/. 

2.3. GALSPECT 

The GALFA-Hl spectrometer, GALSPECT, consists 
of 14 intermediate frequency (IF) to quadrature base- 
band downconverters, 7 dual polarization spectrome- 
ter boards, a CompactPCI computer Board, an Agi- 
lent 8648 A synthesizer, all powered by an uninterrupt- 
ible power supply. The system is designed to take 14 
IF signals split and amplified from the ALFA signal 
chain and generate 14 32-bit spectra each second over 
7679 channels covering ~ 7 MHz as well as 14 32-bit 
spectra each second over 512 channels covering 100 
MHz (see |http://se ti. berkeley.edu/galfa for details). The 
wide-band spectra ar e use d to calibrate the narrow-band 
'science' spectra (see ^4.1| ). The spectrometer has superb 
channel-to-channel isolation due to the use of a Fourier- 
transform polyphase filterbank and loses less than 1% of 
data from 'dead time' during data readout. The spectra 
are tagged with system parameters (telescope pointing, 
IF parameters, etc.) through the AO Shared Common 
Random Access Memory Network (SCRAMNet) and are 
written each second over the AO network to a large server 
for data storage and reduction. 

3. OBSERVING TECHNIQUES 

A suite of observing techniques is used in the GALFA- 
Hl survey, which we describe here. Many of the single- 
beam techniques these mod es were based upon are d e- 



scribed in^'Neil (2002) and Stanimirovic et al. (|2002). 



3.1. GALFA-Hl Survey Philosophy 

When the GALFA-Hl project was started, the impor- 
tance of both the legacy data product and target-oriented 
observing approach were recognized. These two goals 
could easily come into conflict if not treated carefully: 
too much freedom for observers to pursue individual 
projects could yield data that would be difficult to com- 
bine into a final survey product, too many constraints on 
the observations would make scientific goals hard to pur- 
sue. To facilitate the creation of a final survey, all of the 
GALFA-Hl data are taken in a single data-taking mode, 
with a fixed bandwidth, resolution and frequency setting. 
This is made possible by the GALSPECT spectrometer, 
which simultaneously has narrow enough channels (0.18 
km s~^ ) for the most stringent Galacti c Hi science crite- 
ria (Hi narrow-line self- absorption, e.g., Krco et al.|2007[ ) 
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survey 


telescope 


era 


area 


resolution 


velocity resolution 


sensitivity 


LAB 


Dwingeloo and lAR 


1989 - 2005 


100% 


36' 


1.3 km s-i 


-^100 mK 


IGPS 


various 


1995 - 2010 


5% 


1^-2' 


1.3 km s-i 


2700 mK - 1500 mK 


GALFA-Hi 


Arecibo 


2005 - 


32% 


4.0' 


0.2 km s-i 


60 mK - 140 mK 


GASS 


Parkes 


2005 - 2009 


50% 


15' 


1.0 km s-i 


60 mK 


EBHIS 


Effelsberg 


2009 - 


~ 50% 


9' 


1.4 km s"-*^ 


90 mK 


GASKAP 


ASKAP 


2012 - 


25% 


20" - 3' 


0.2 km s-i 


3000 mK - 40 mK 



TABLE 1 

A COMPARISON OF THE MODERN GALACTIC Hi SURVEYS: LAB JKALBERLA ET AL."2005'), IGPS, ('TAYLOR ET AL."2003"StIL ET AL.| |2006| 

AND McClure-Griffiths ET AL. 2005), GALFA-Hi, GASS (McClure- Griffiths et al. 2009), EBHiy (Winkel et al. 2U1(^ and 
GASKAP. Era describes the time period over which the data were taken and reduced. Area is shown in percent of entire 

SKY COVERED. THE IGPS DATA WERE TAKEN AT DIFFERENT RESOLUTIONS AND SENSITIVITIES BY DIFFERENT SUB-SURVEYS. THE GASKAP 
DATA CAN BE REDUCED IN DIFFERENT WAYS, LEADING TO DIFFERENT RESOLUTIONS AND SENSITIVITIES. THE SENSITIVITY LISTED IS FOR A 1 

KM S~^ CHANNEL. 
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Fig. 1. — The integration time of data release 1. The top panel is in equatorial coordinates, the bottom panel is in Galactic coordinates, 
with integration time per beam in grayscale. The rectangle in the top image indicates the total skyavailable to the Arecibo observatory, 
as does the black line in the bottom image. The dashed line represents the area shown in Figure [2] the dash-dotted line represents the 
sub-area shown in Figure [To] 



Fig. 3. — A model of the response function of each of the 7 
ALFA beams. Note the asymmetric sidelobe response in beams 1 
- 6. The beams are separated from each other in this diagram more 
than they are on the sky to show the sidelobes without coUision. 



as well as enough bandwidth (+700 to -700 km s~^ ) 
to easily encompass the fastest know n moving Galacti c 
Hi gas (very high velocity clouds, e.g., Peek et al.||2007 ). 
In this way all projects can be served by the same spec- 
trometer setting. It is also important to allow individual 
observers to image higher-priority projects first. To ob- 
serve the totality of the Arecibo sky, even at the shallow- 
est possible level, would require more than 1700 hours of 
observing, thus delaying the publication of particular sci- 
ence of interest. The decision to allow individual Pis to 
choose science targets has indeed facilitated the timely 
release of data and the inclusion of students in the^ survey 
process (Goldsmith & GALFA-Consortium 2003). The 
DRl coverage maps (Figure[T]) show the effect of this phi- 
losophy, with each area having been observed for varying 
lengths of time and in varying observing modes (see Ta- 
ble |2|. As the survey continues, the entire Arecibo sky 
will DC mapped, but at varying levels of sensitivity. 

3.2. Least-Squares Frequency Switching 

In radio heterodyne spectroscopy, the measured spec- 
trum is the product of the radio-frequency (RF) power 
and the intermediate-frequency (IF) gain spectra; to ob- 
tain the RF power spectrum, one must divide the on- 
source measured spectrum (the 'on spectrum') by the IF 
gain spectrum. To accurately determine IF gain, we em- 
ploy a new method call ed least-squares frequency switch- 
ing (LSFS; Heiles"2005b). LSFS is fundamentally differ- 
ent from frequency switching. In frequency switching 



(the methodology us ed by GASS, for example; McClure- 
Griffiths et aT][2QQ9l ), one takes two spectra — one on-line 



and one ott-line spectrum — and (1) finds the difference, 
to remove the continuum (frequency-independent) sys- 
tem temperature, which is large; and (2) divides the dif- 
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data set 


principal investigator(s) 


time 


area 


start date 


end date 






hours 


sq. deg. 


MM/DD/YY 


MM/DD/YY 


DRl 


3046 


7520 


05/16/05 


02/26/09 


DRl-F 


1429 


4043 


05/16/05 


02/26/09 


a2032 


Stanimirovic 


41 


203 


05/16/05 


09/03/05 


a2011 


Peek &: Heiles 


13 


111 


05/21/05 


06/12/05 


a2050 


Peek &; Heiles 


41 


414 


06/25/05 


01/21/06 


a2004 


Paul Goldsmith 


38 


337 


07/09/05 


12/31/07 


TOGS 


Putman &; Stanimirovic 


930 


1766 


08/17/05 


01/14/09 


a2172 


Heiles &; Peek 


106 


974 


06/07/06 


02/18/07 


a2222 


Peek 


63 


615 


07/29/06 


01/17/07 


a2174 


Lewis Knee 


33 


316 


09/07/06 


11/05/06 


TOGS2 


Stanimirovic 


164 


1526 


11/14/08 


02/26/09 


DRl-S 


1616 


3753 


06/17/05 


07/04/08 


a2011 


Peek 


28 


285 


06/17/05 


10/22/05 


TOGS 


Putman &; Stanimirovic 


1367 


2672 


12/10/05 


05/31/08 


a2221 


Leonidas Dedes 


10 


104 


08/19/06 


09/17/06 


a2220 


Heiles & Peek 


40 


429 


10/28/06 


12/04/06 


a2060 


Putman 


171 


864 


11/03/06 


07/04/08 



TABLE 2 

A TABLE OF PARAMETERS FOR DRl AND ITS SUBSETS. DRl IS COMPOSED OF TWO INDEPENDENTLY REDUCED SUBREGIONS, DRl-S AND 
DRl-F, ALL SHOWN IN BOLD. EACH OF THESE SUBREGIONS ARE THEMSELVES COMPOSED OF SUBREGIONS, REPRESENTING INDIVIDUAL PI 
PROJECTS THAT HAVE BEEN COMBINED TOGETHER. THEY ARE INDICATED HERE BY THEIR ArECIBO PROJECT NUMBER. THE "TURN ON 

GALFA Survey", a fully commensal set of observations that make up the lion's share of DRl, is indicated by TOGS 
(commensal with ALFALFA) and TOGS2 (commensal with GALFACTS). 



ference by the off-line spectrum, which corrects for the 
frequency-dependent system gain. This technique intro- 
duces two artifacts. First, the system temperature is not 
independent of frequency, and in particular at Arecibo 
it has variations on the ^ 1 M Hz scale, known as fixed- 
pattern noise (FPN; see ^4.3 and Figure |4|. Thus, when 
taking the difference in (Ij, the FPN in both the on- 
line and the off-line spectra appears in the difference 
spectrum; these are uncorrelated, so the FPN amplitude 
increases by ~ V2. Second, the system gain depends 
mainly on IF frequency, and the division in (2) reliably 
corrects for this IF gain component. However, the sys- 
tem gain also varies with RF frequency to some extent, 
and (as with the FPN) RF gain variations in the off-line 
spectrum appear in the final result. Our LSFS technique 
avoids some of these problems because it avoids the need 
for subtracting an off-line spectrum, and it derives the IF 
gain alone, so that the final derived spectral gain errors 
come only from those in the RF spectrum. 

In LSFS, one sets the local oscillator (LO) frequency 
to a number of different values so that the RF power 
spectrum and the IF gain spectrum can be evaluated as 
distinct entities using a least-squares technique ( Heiles, 
|2007). In practice, at the beginning of each observation 
we spend 10 minutes recording data at a single sky po- 
sition at a number of different frequencies. We used the 
minimum-redundancy schema with 7 LO settings. In this 
"MR7" schema, the LO frequency separations range from 
1 to 39 times the minimum separation A/. One needs the 
maximum separation to be a reasonable fraction of the 
total spectrometer bandwidth so that the full-bandwidth 
gain spectrum is accurately defined. With our large num- 
ber of channels, this means that A/ must be much larger 
than the channel spacing Sf. We chose = 225, which 

gave a maximum separation of 90% of the full spectrom- 
eter bandwidth. 
When > 1, one needs a reliable way to degrade 



the resolution of the original 7679-channel spectra to 
reduce the number of channels in the LSFS technique 
to a manageable number. One applies the LSFS tech- 
nique to these degraded-resolution spectra to derive a 
degraded-resolution IF gain spectrum. Finally, one in- 
terpolates the degraded-resolution gain spectrum to re- 
gain the original number of spectral channels. We used 
standard Fourier series to both degrade the original res- 
olution and to then recover it, reducing the number of 
channels in the degraded-resolution spectra to 512. 

3.3. Scanning Methods 
3.3.1. Drift Scanning 

The simplest method of observing is the drift-scanning 
method. Indeed, this method traces back to the dawn 
of radio astronomy as the telescope is simply left at a 
fixed position at transit (in our case at azimuth = 180 or 
azimuth = 0) and the sky is allowed to drift by (Reber 
1940). To space the beams equally 2') the AEFa 
beams are rotated ~ 19° from the orientation shown 
in Figure [3j as shown in Figure |5j This scanning pat- 
tern holds the maximum number of telescope parame- 
ters fixed, which minimizes variability in baseline ripple. 
This pattern also has a number of disadvantages. Firstly, 
it is unwieldy for regions that are relatively large in Dec- 
lination (or small in Right Ascension). Secondly, since 
the scans are all parallel on the sky, there is no easy 
way t o cros s-calibrate the scans, as they do not intersect 
(see ^3.3.2). Thirdly, the telescope has limitations on 
the zenith angle at which it can point, so regions near 
zenith are inaccessible in the standard mode. It is possi- 
ble to drift scan in positions off transit to get Declination 
range near zenith, but one sacrifices constancy of FPN. 
This method is used solel y by the Areci bo Legacy Fast 
ALFA Survey (ALFALFA; ICiovanelh et' al. 2005), which 
is a wide-area survey with which GALFA-Hl commen- 
sally observes. Data taken commensally with ALFALFA 




1418.9 1421.9 1418.9 1421.9 1418.9 1421.! 

Frequency [MHz] 

Fig. 4. — An example of the FPN. The left panel shows the uncor- 
rected FPN in each of the seven beams, averaged over polarization 
and over a few hours of observation. A gray line at zero is shown 
for comparison. The middle panel shows these same spectra less 
their average. Plotted behind these average spectra is the fit to 
these spectra in gray, as in the ripple component of Equation ^ 
The right panel shows the original spectra less the fit in the center 
panel. The ripple has been significantly reduced. 



(the Turn On GALFA Survey, or TOGS, project) repre- 
sent 75% of the raw data that make up DRl, and cover 
60% of the DRl sky. This means that the TOGS areas 
are some of the highest sensitivity areas of the DRl data 
set (see the dark bands at fixed dechnation in Figure [T]). 

3.3.2. Basketweave Scanning 

The preferred mode for GALFA-Hl observations is bas- 
ketweave or 'meridian-nodding' scanning. In this mode 
the telescope is kept at the meridian and is driven up and 
down in zenith angle over the chosen Declination range. 
Each day the starting position of these nods is displaced 
in Right Ascension by ~ 12' such that the entire re- 
gion is covered, as in Figure [5] Alternatively, the scan 
can be started at a fixed RA, and the Declination can 
be varied, which is t he method used by our commensal 
partner GALFACTS (|Guram & Taylor|2009|) . There are 



six distinct drive speeds (or 'gears') that generate equally 
spaced tracks with each of the seven beams: three with 
an ALFA rotation angle of 0° and three with an ALFA 
rotation angle of 30°. Practically, only the two fastest 
gears, one with rotation angle 0° and one rotation angle 
of 30°, are used. This observing method has the ad- 
vantage that it generates a large number of intersection 
points among beam tracks at which the relative gain of 



QJ 
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Basketweave 




Drift 



Right Ascension 



Fig. 5. — The two modes of scanning used by GALFA-Hi Data 
Release 1. In the examples, the empty circles represent the initial 
position of ALFA for a single day's observation and the filled circles 
represent the final position of ALFA for the day's observation. The 
dashed lines represent complementary scans taken on other days 
of observations. 



the beams can be measured (as in Haslam et al. 1981). 
This allows us to greatly reduce the striation in the map 
caused by gain variation. It does this while maintaining a 
fixed azimuth, which significantly limits the fiuct nations 
in FPN. 

A somewhat modified version of the basketweave can 
be employed at Declinations that pass near zenith at 
Arecibo where the traditional basketweave scan cannot 
be used. In this method, called 'azimuth wagging', the 
same scan pattern on the sky is achieved, but by posi- 
tioning the telescope at near 90 or 270 degrees azimuth. 
This method is inferior to the standard basketweave in 
that the azimuth, zenith angle and ALFA rotation angle 
must all change continuously during the observation. 

4. THE REDUCTION PROCESS 

4.1. Bandpass Fitting with LSFS 

As discussed in §3.2[ we need to correct each spectrum 
for the IF gain. At the beginning of each observing pe- 
riod, we run an LSFS calibration to obtain this IF gain 
spectrum and we use this spectrum for the whole period's 
data. This works because the IF gain spectrum is largely 
stable. We find that from day to day the IF solutions typ- 
ically vary with amplitudes of only dozens of mK, but 
larger jumps of up a few K are observed occasionally. 
We attempt to maximize this stability by setting all as- 
pects of the IF electronics, including all attenuator and 
amplifier gain values, to standard values. Occasionally 
LSFS data are contaminated, either by RFI or by im- 
proper implementation during the observation. In this 
case we use the LSFS data product from the previous or 
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following day of observation. 

4.2. Initial Calibrations 

After the effects of the bandpass are removed, the data 
are searched for possible contamination by reflections in 
the signal chain (see Figure |6|. Any impedance mis- 
match between cables or signal transmitting hardware 
can cause echoes within the system. Similar to the re- 
flections in the superstructure and geodetic dome of the 
telescope, these reflections in the cables produce excess 
power at lags equivalent to the reflection's differential 
travel time, and therefore cause ripples in the baseline of 
the spectrum. These reflections have the very nice qual- 
ity that they are each only a single Fourier mode: as the 
reflection geometry is only one dimensional, there is no 
room for complex reflection patterns as there is in the 
telescope superstructure. This ripple is therefore very 
easily removed by examining the data in Fourier space 
and simply removing the clear reflection tone. We also 
remove any single-channel RFI that is detectable above 
the noise in the data smoothed over a ten-second inter- 
val by setting it to the median of the adjacent channels. 
A secondary RFI effect is also searched for: radar in- 
termodulation products. In certain LO configurations, a 
periodic radar signal emanating from the San Juan air- 
port can generate periodic pulses in our data that can be 
hundreds of channels wide and as much as 10 K in am- 
plitude. These pulses have a very well known period and 
spectral shape, so we can search for them and eliminate 
them from the data relatively easily. After these spectral 
cleaning procedures are done the data are shifted into the 
local standard of rest (LSR) frame. 

4.3. Fixed-Pattern Noise Subtraction 

It was known since long before the GALFA-Hl project 
was undertaken that spectra taken at the Arecibo 305 
m (and all radio telescopes) are typically cont aminated 
by fix ed pattern noise (FPN; see, for example, Padman" 
[1977^ and references therein). FPN (also called 'baseline 
ripple') at Arecibo manifests as a low- level wave that per- 
meates the entire spectrum, with peak-to-trough varia- 
tions of ^ 400 mK across a single, polarization-averaged 
band (see the first panel in Figure |4l for a detailed de- 
scription of ALFA's FPN see Heiles'T2QQ5a])). Note that 
while this is certainly strong FPN, it does not dramati- 
cally exceed values found for telescopes with completely 
unblocked apertures; the Green Ban k Telescope has FPN 
at about 30% this level in L-band (Fisher et al. 2003). 
At Arecibo, FPN contamination to the spectra is caused 
by reflections inside the Arecibo 305 m superstructure 
and geodetic dome. Any signal that takes more than a 
single path to the receivers will cause a spurious increase 
in the amplitude of the auto-correlation function (ACF) 
in the channel corresponding to the delay between the 
paths. In the case of a signal that takes both a stan- 
dard path to the receivers and one that has an extra 
reflection between the geodetic dome and the main re- 
flector, for instance, the path length difference is on the 
order of 300 meters, which corresponds to a delay of 1 /is. 
This 1 /is delay causes an increased ACF at the 1 /is lag 
which, in turn, generates a ripple with a period across 
the spectrum of 1 MHz, or 200 km s~^. Since Galac- 
tic Hi typically covers a comparable velocity range, this 



FPN is capable of corrupting the final data product if 
not addressed adequately. 

Removing the FPN caused by the reflections in the 
telescope superstructure and geodetic dome is much more 
difficult than removing the effects of reflections within 
the signal chain. This is because the FPN is not dom- 
inated by a single Fourier component, but rather by a 
large range of Fourier modes, with associated delay rang- 
ing from 0.5 to 2 /is. The FPN was examined in some 
detail in 'Heiles (2005a). This work showed that the rip- 
ples are 100% polarized, uncorrelated from receiver to 
receiver, relatively static over time and dependent upon 
ALFA rotation angle, telescope azimuth and zenith an- 
gle. 

To reduce the effects of the FPN we use the fact that 
each beam has an uncorrelated and relatively fixed pat- 
tern over the course of an observation in our two main 
observing modes, drift scanning and basketweave scan- 
ning. In the case of drift scanning all telescope position 
parameters are held fixed, so we expect no variation in 
the FPN. In the case of basketweave scanning only the 
zenith angle is changed, which causes an average change 
of ~ 25% in the FPN over the full range of zenith angle, 
which we ignore. After averaging the two polarizations 
within a beam and applying a simple overall gain calibra- 
tion, we reduce the effect of the FPN by first generating 
the average spectrum, Tavg {^) for all 7 beams over the 
course of a single day's observation. We then fit the dif- 
ference between Tavg (^), and the average over beam n, 
Tn{iy), with 

m 

Tn (ly) - Tavg (^) = ^ {aiSifl (tiU) + biCOS (tiiy)) (1) 



^1 ASnd{u-Uj) 



where, ti ranges over the observed FPN frequencies (0.5- 
2 /is) and Uj ranges over the frequencies where Galactic 
disk Hi is typically observed, d is equivalent to the Dirac 
delta function, to avoid confusion with Declination, 6. 
The first component of the fit is the baseline ripple, the 
second component is the residual Hi contribution. By 
fitting the residual Hi contribution with this Taylor ex- 
pansion in a and S we are able to accurately capture the 
Hi residuals while only fitting 2 parameters per chan- 
nel, rather than the 6 required if we were to fit each 
beam individually. By subtracting only the baseline rip- 
ple component of the fit for each beam, we are able to 
reduce the amplitude of the baseline ripple in each beam 
by ~ a/7 without disturbing the Hi data. An example of 
this method is shown in Figure |4] 

4.4. Intersection- Point Calibration 

To make a correct map from data taken with a single 
dish radio telescope, the overall gains need to be accu- 
rately calibrated. These gains can vary significantly from 
day to day, based on both the performance of the tele- 
scope and the receiver (see Figure [7| which shows the 
variation on day and minute timescales). These gains 
may even vary during long observations. It is possible to 
calibrate these variable gains with the calibration noise 
diodes in the ALFA system, but the amplitude at the 
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Fig. 6. — The plots on the left are time averaged spectra from a drift scan observation in the original state (top) and with the cable 
reflection removed (bottom). The plots to the right are the Fourier transform of the spectra, with the region where the reflection was 
removed highlighted. Note that the digital Fourier transform of a monochrome signal in general has both a central spike and "ringing" 
component in Fourier space; our fltting also largely reduces the ringing component. The longer period and more chaotic FPN has not been 
reduced in this example. 



receiver that these diodes generate can also vary signifi- 
cantly from day to day. Uncalibrated gains lead to stri- 
ation in the maps and inaccurate Hi amplitudes. 

To determine the relative gains of each beam over each 
day of observation, we use the constraint that at places 
where beam tracks intersect they must observe the same 
region of sky, so any differences in their spectra must be 
dominated by variation in the gain parameters. This 
method is most effective in basketweave scans where 
there are many intersection points - a region 2 hours 
across in Right Ascension and 17 degrees in Declination 
has almost 10^ intersection points. If we wish to deter- 
mine the relative gains in a set of drift scans, we must 
conduct other observations that intersect the drift scans, 
as drift scans are parallel on the sky. 

Once the intersection points are found, we determine 
the relative gains at these points by comparing the am- 
plitudes of the observed Hi spectrum in each beam. Each 
beam is assumed to have an overall gain for each day's 
observation plus some variability with time during the 
observation that we parameterize as a sum of Fourier 
components. Since only relative gains are measured, the 
overall gain solution must average to unity. We construct 
a least-squares fit to the measured relative gains at in- 
tersection points for the various day and time dependent 
gain parameters and then apply these parameters to the 
entire data set to remove the relative gains. 



4.5. Gridding 

Once the spectra, or time-ordered data (TOD), have 
been fully calibrated, we apply them to a grid in the 
image plane. This grid can be made in any proje ction 
in the World Coor dinate System as described by Cal-| 
abretta & Greisen'' (2002), but we project to Cartesian 
coordinates in Right Ascension and Declination with an 
orthogonal graticule. We use a Gaussian sampling ker- 
nel with FWHM = 3' to place flux on the grid; this size 
sampling kernel smoothes out any pixels that may have 
been missed by the observing pattern without degrading 
the resolution too significantly. Our method is anchored 
by the IDL code SDGRID.PRO, a code based upon the 
IDL gridzilla implementation by Snezana Stani mirovic. 
SDGR ID.PRO uses the core logic described in |Dickey| 
( 2002[) a nd many of the met hods developed for Parkes 
gridzilla (Barnes et al. 2001) and the AIPS procedure 
SDGRID. At th is pomt the first sidelobe correction is ap- 
plied (see ^4.6, below). Once sidelobe correction is com- 
plete, a final amplitude correction is applied to make the 
observations consistent with the lower resolut ion Leiden- 
Argentine-Bonn Survey ( Kalberla et al.|20Q5 ). We select 
a large region of sky towards the plane, where stray radi- 
ation will be relatively weak, and compare the GALFA- 
Hl spectrum to the LAB spectrum for the region. This 
comparison results in a 10% to 15% adjustment of the 
final GALFA-Hl data cubes. 
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suming the sky is identical to our data cube from ^4.5 



20 30 
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Fig. 7. — The relative gains between two consecutive days of ob- 
servation, taken at identical positions. Each coincident integration 
is compared for each polarization-averaged beam, plotted in black 
(one line for each beam). We find the relative beam gain between 
the two days at each second by fitting to the Hi spectrum taken on 
the first day plotted against the Hi spectrum taken on the second 
day, where the relative gain is the slope of the line. Each relative 
beam gain is averaged over 100 seconds to smooth the gain trends. 
The smoothed beam gain data are plotted in black, with the un- 
smoothed beam relative gains in light gray for comparison. The 
noise in the unsmoothed gain fits changes over the course of the 
observation time as the observations leave the Galactic disk and 
the Hi amplitude decreases significantly. 



4.6. First Sidelobe Calibration 

Each ALFA beam has a significantly different beam 
shape, because of the very fast optics of the telescope 



( Heiles][2Q Q4 ) . The center beam (beam 0) has a well be- 
haved beam shape, with sidelobe efficiency, r]sLj near 
10%. Each other beam (beams 1-6) has significant first 
and second sidelobes which depend upon the relative dis- 
placement of the beam to beam (see Figure [3|. The 
first sidelobes are ^ 5' offset from the main beams with 
r]SL near 20% for beams 1 to 6. Since we do not Nyquist 
sample the sky with each beam, but rather with a com- 
bination of all the beams, we end up with a map that has 
an ill-defined associated beam pattern. This means that 
many of the standard methods devised for removing the 
effect of beam shape fro m maps, such as the maximum 
entropy method (Cornwell & Braun||1989|and references 
therein) and CLEAN (Hogbom 1974), cannot be used 
and we must devise our own method for removing the 
effects of these sidelobes. We note that the GALFACTS 
survey has de veloped a 'multi-beam CLEAN' algorithm 
(Guram & Taylor 2009), though it is primarily designed 
tor point-source reconstruction, rather than the diffuse 
Galactic Hi signal we detect. 

To counteract the effect of these inconsistent, asym- 
metrical beams we use a strategy similar to that em- 
ployed to correct the Leiden-D win^aloo Survey, as out- 
lined in Hartmann et al.'('1996). Because we have seven 
different beams, rather than a single convolving beam, 
we must modify their method somewhat. The method, in 
essence, is to determine a "contamination spectrum" for 
each of the TOD; these contamination spectra represent 
what we expect the sidelobes alone would detect. We 
construct these contamination spectra by 're-observing' 



These contamination spectra are then subtracted from 
the original spectra, 

T'{v,t)=T{v, t) /r]MB - VslTc [y, t) /tjub , (2) 

Where T {v^ t) are the original TOD, Tq (i^, t) are the 
contamination spectra, and 77MB and 77SL are the main 
beam and sidelobe efficiencies, respectively. These cor- 
rected spectra, T' {v^ t) are then re-gridded to make a fi- 
nal map. While there are second-order effects that stem 
from using the contaminated map as the true sky im- 
age, we have found that these effects are small, of order 
T^gL, and typically smaller than the effect from unmod- 
eled higher order sidelobes. 

We measure the effective beamsize of the final maps by 
measuring the size of compact radio sources in Hi absorp- 
tion. We examine the 100 brightest NVSS sources that 
are both within DRl and in regions of the sky where 
the Galactic Hi does not vary significantly over the re- 
gion near the source. We find that an effective beam- 
size of 4.1' in Declination and 3.9' in Right Ascension. 
This ellipsoidal nature of the effective beam is due to the 
known elongati on of the Ar ecibo beam along the zenith 
angle axis (e. g. Heiles 2004), which is kept oriented along 
Declination durmg our observations. 

5. DATA PROCESSING: AN EXAMPLE REGION 

As an example of the GALFA-Hl survey observations 
and data processing methods we have described in the 
previous sections, we show here a specific region ob- 
served over several months in 2006. The entire region 
is 20° X 17°, centered on {a, 3) = (17^, +8°), and was 
chosen by the investigators (Putman and Stanimirovic) 
because HVC complex C is visible in the region. It was 
observed with the basketweave method (see [ 3.3.2 ) twice, 
to achieve a higher signal-to-noise. 

In this example we investigate a single survey cube, 
which is 512 x 512 V square pixels and in the standard 
Cartesian projection and 2048 frequency channels. Fig- 
ure [8] shows the average spectrum over the cube and 
the two velocity ranges we consider. We isolate these 
velocity ranges to demonstrate the effect of each stage 
of the data reduction on both strong Hi (on-line) and 
weak Hi (off-line) maps. We demonstrate the effects of 
each stage of reduction in Figure [9) by showing a map 
of the uncalibrated data at the top left, fully-calibrated 
data at the bottom right, and the corrections applied at 
each step in between. We set these up like an equation, 
with the final data cube equalling the original data cube 
less each isolated contaminant ("initial calibrations and 
FPN" , "crossing point gain calibrations" , and "first side- 
lobe corrections"). The maps are scaled to fit the range 
of observed surface brightnesses, from K to 2 K in the 
case of off-line map (top) and from 7 K to 23 K in the 
case of the on-line map (bottom). Correction maps are 
scaled from -0.1 to 0.1 K (top 3 correction maps) and 
-1 K to 1 K (bottom 3 correction maps). Note that the 
software does not typically generate these intermediate 
images, but we show them here for clarity. 

The first two panels ("original, uncorrected image") 
in the top and bottom section of Figure [9] show the 
integrated intensity map of the off-line and on-line re- 
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Fig. 8. — The average spectrum over the region discussed in ^ 
The on-hne and off-hne maps were integrated over the velocity 
ranges indicated in light gray and dark gray domains, respectively. 



gions without any cali bration past the initial calibra- 
tions described in §4.2[ The strong striation along the 
scan patterns from poorly calibrated gains and baseline 
ripple are quite evident. The following three pairs of 
plots show the initial gain calibrations and FPNcorrec- 
tion (^4.3), crossing-point gain calibration (^4.4), and 
the first-sidelobe correction ( ^4.6[ ), respectively. For the 
first (FPN) correction, there is a strong effect in both 
the on- and off-line data. In the case of the off-line plot 
(top), the cleaning is primarily due to the ripple reduc- 
tion itself, as the amplitude of the features on the sky 
is similar to the amplitude of the ripple. In the on-line 
slice (bottom) the cleaning is not due to baseline ripple, 
but rather the fact that a rough gain calibration applied 
when the baseline ripple is corrected. In the second pair, 
"crossing point gain calibration" , there is only a modest 
change due to gains, whereas the much brighter on-line 
data show significant changes in some regions. The fi- 
nal step of first-sidelobe correction has most effect at 
the bright peaks of the slices where data had been er- 
roneously scattered, and has a noticeable effect on both 
off-line and on-line slices. The last pair of maps shows the 
same region and velocity range as the first set of maps, 
but with all three data reduction methods applied. The 
off-line data are significantly improved, but given the low 
dynamic range of the map (~ 1 K), it is not surprising 
that residual ripple on the order of 0.1 K (see Figure 
|4| is still visible. The on-line data are also drastically 
improved, and even with very close inspection it is very 
hard to determine the presence of any systematic arti- 
facts. These maps demonstrate the very high quality of 
the GALFA-Hl data products. 



6. DATA RELEASE 1 

Data Release 1 (DRl) is the first full data release of the 
GALFA-Hl survey. The 3046 hours of data from 12 inde- 
pendent projects were taken from May 2005 to February 
2009: details can be found in Table [2] Integration times 
per beam area range from ~10 seconds to ~60 seconds 
(see Figure [T]) , as different regions were observed in dif- 



ferent modes, and sometimes covered by more than one 
project. These integration times correspond to 140 mK 
and 60 mK RMS noise, respectively, when smoothed to 
a 1 km s~^ channel width. The median integration time 
per beam is 30 seconds (80 mK). DRl was reduced in 
two separate stages; DRl-S (spring), covering the sky 
towards 12h RA, and DRl-F (fall), covering the sky to- 
wards Oh RA. DRl-S was made publicly available in Jan- 
uary 2009, and the entirety of DRl is released with this 
work. The two areas have an overlap of 180 square de- 
grees toward {a^S) = (6^,25°), and give consistent re- 
sults. In Figure [lOj we demonstrate the quality of the 
DRl data set, and the fantastic structure of the ISM at 
4' resolution, by showing three sections of velocity space 
over the same arbitrary chosen 40° x 18° 10' region. The 
intense filamentary and diversity of cloud structures is 
immediately apparent. 

6.1. DRl data reduction and known systematics 

Data reduction was very similar for both halves of 
DRl, but some differences exist. For DRl a routine for 
removing radio intermodulation products was added for 
use with T0GS2 (commensal with GALFACTS; see Ta- 
ble |2|, because of the specific LO setup used. The exact 
protocol for data flagging was improved between DRl-S 
and DRl-F, which has led to improved data quality in 
DRl-F. Primarily residual fixed pattern noise is lower in 
DRl-F than in DRl-S, owing to improved algorithms and 
more precise by-eye checking procedures. An example of 
the scale of strong FPN residuals in the two halves of the 
data set can be seen in Figuref 



11 There are also fewer oc- 



casional RFI bursts and glitches in DRl-F, though these 
are rare in either data set. 

The Arecibo 305 m is also known to have significant 
distant sidelobes, or stray radiation contamination. This 
effect can be important in Galactic Hi work, as the sky 
is bright at Galactic velocities in all directions. Since 
the stray sidelobes are large and low-gain, the effect is 
typically to add a wide, low amplitude signal that varies 
slowly with position. We intend to implement some stray 
radiation correction with future data releases, but for 
DRl we simply attempt to parameterize the amplitude 
of the systematic error. Figure [l2| shows spectra towards 
two positions in the sky, one oFthe lowest column parts 
of the DRl sky (/ = 210°, b = 65°), and one towards a 
low column part of the DRl sky with the largest amount 
of bright Galactic Hi high above the horizon, to maxi- 
mize stray contamination (/ = 45°, b = 45°). We show 
both the average LAB spectrum for the region and the 
difference between the GALFA-Hl average spectrum for 
the region and the LAB spectrum for the region to high- 
light the effect of stray radiation. We find that the stray 
effect is strongest closest to the line (|Flsr| < 50 km 
s~^), with an amplitude below 200 mK outside ± 20 km 
s~^ LSR, and that it is strongest when the Galactic plane 
is high in the sky. Within ± 20 km s~^ LSR, we cannot 
accurately enough match the beam shape and sensitivity 
profile of the LAB data to get a reliable measurement of 
the stray radiation effect, though it is unlikely to exceed 
500 mK. 

7. FUTURE DIRECTIONS AND UNRESOLVED 
ISSUES 
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Fig. 9. — An example region of the sky towards (a, 6) = (17^ 52™, +10° 21^). Tlie top and bottom areas represent the dark (off-hne) and 
hght (on-hne) regions in Figure [s] respectively. The maps are scaled to fit the range of observed surface brightnesses, from K to 2 K in 
the case of off-line map (top) ana from 7 K to 23 K in the case of the on-line map (bottom). The calibration images are scaled from -0.1 
K to 0.1 K in the top panels, and from -1 K to 1 K on the bottom panels. See JsTfor details. 
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Fig. 11. — Panels representing typical-to-strong residual FPN 
in the DRl-S (top) and DRl-F (bottom) data. The left panels 
are data cubes averaged over a few degrees in RA in high-latitude 
regions with typical strong FPN taken in drift mode. The data have 
had variation above a scale of 15^ removed to highlight the beam- 
to-beam FPN variation. The right panels show the 1-cr (lower, 
smooth) and peak-to-trough (higher, jagged) amplitude of the FPN 
as a function of velocity. The data show that DRl-S is somewhat 
worse than DRl-F, that the FPN errors increase somewhat as the 
velocity nears the Hi line, and that the typical FPN amplitude is 
150 mK peak-to-trough, or 20 mK 1-cr in these more contaminated 
regions. 



We have made significant strides towards reliable, 
artifact-free data, and we expect that these data will be a 
useful resource for the astronomy community. To further 
this goal we intend to continue to pursue known sources 
of systematic error in our data, including the effects of 
more distant sidelobes and stray radiation. 

7.1. Further Sidelobe Calibration 

While the primary source of stray radiation is un- 
doubtably the first sidelobes, the ALFA beam pattern 
is known to have significant power in secondary and 
higher-order sidelobes. These sidelobes roughly follow an 
Airy pattern m odel with elliptical cross section (Cortes^ 
Medellm||2QQ2 ). Second sidelobes are more distant from 
the center of tne observation, so while the contamination 
amplitude will be lower, the signal in the second sidelobe 
will be typically be more divergent from the main beam 
than in the first sidelobe. In cases where there is very fast 
variation in the sky signal this second sidelobe can have 
a significant effect, and we believe that it is currently the 
dominant contribution to uncorrected sidelobe radiation. 

To correct for the second and higher-order sidelobes 
these sidelobes must be carefully measured. Some raster 
maps have been made of the second sidelobe response, 
but they are significantly incomplete. To parameterize 
the basic components of second sidelobe w e will measur e 
it with the 'spider scan' method outlined in Heiles (2004). 



We will then add the second sidelobe response to our 
calibration method described in §4.6[ 

7.2. Stray Radiation Calibration 

We have shown that there is non-trivial stray radiation 
in our data. While the stray radiation does not pose a 
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Fig. 12. — Spectra toward {I, b) = (210°, 65°), a very low column 
region in DRl (bottom) and (l,b) = (45°, 45°), a low column region 
observed while the bright Galactic plane is well above the horizon 
(top). The thick gray line is the LAB data, and the thinner black 
line is the GALFA-Hi data minus the LAB data, which outside 
regions of bright Hi represents the contamination by stray radiation 
towards these directions in the GALFA-Hl data. We have marked 
with a gray box the velocity range where we expect the difference 
to be dominated by beam mismatch between the LAB survey and 
our reconstructed LAB beam with GALFA-Hl data, and where we 
can therefore not use the plotted difference as an indicator of stray 
radiation strength. 



problem to conduct most of the science we are pursuing, 
including any science dependent upon small scale mor- 
phological features, or features at high velocity, it can 
be a problematic contaminant, especially when precise 
total Hi columns are needed. We plan to correct for this 
stray effect in future data releases by bootstrapping our 
results to the LAB data, and removing the differences we 
detect. This will allow us to treat the stray data prop- 
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erly without modeling the very complex feed response far 
from the beam center generated by the blocked Arecibo 
305 m aperture. 

8. CONCLUSION 

In this work we described the systems observ- 
ing modes and data reduction methods (M| used to 
produce the CjALFA-Hi survey. We demonstrated the 
efficacy of these new and inventive techniques as well as 
their current limitations, and showed the high quality 
of the resulting data products (^. We presented the 
GALFA-Hl Data Release 1 and discussed methods 
for improving our final data product in the future (Q. 
We found that we were capable of dramatically reduc- 
ing the impact of these systematic effects, in many cases 
making them undetectable. 

The DRl data product can be found at 
https://purcell.ssl.berkeley.edu/, downloadable either in 
the original survey data cubes or in any user-defined 
region and spectral resolution. The original data cubes 
are 512' in small-circle Right Ascension and Declination, 
with 2048 channels, either at the original resolution of 
0.184 km s~^ with a velocity range from —188 km s~^ to 
188 km s~^ LSR or at a degraded resolution of 0.736 km 
s~^ over the full bandwidth. 

GALFA-Hl is very much an ongoing survey with hours 
of data still being taken daily through commensal and 
dedicated projects. As such, future data releases are 
planned including these new data and improved reduc- 
tion methodologies, and data release 2 is tentatively 
planned for late 2011. Many astronomers around the 
world are engaged in active research with GALFA- 



Hl data, and we plan to continue to provide high quality 
data to these astronomers as well as the rest of the as- 
tronomical community in the future. 
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